Ultracold atomic gases in non-Abelian gauge potentials: The case of constant Wilson loop 
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Nowadays it is experimentally feasible to create artificial, and in particular, non-Abelian gauge potentials for 
ultracold atoms trapped in optical lattices. Motivated by this fact, we investigate the fundamental properties of 
an ultracold Fermi gas in a non-Abelian U (2) gauge potential characterized by a constant Wilson loop. Under 
this specific condition, the energy spectrum exhibits a robust band structure with large gaps and reveals a new 
fractal figure. The transverse conductivity is related to topological invariants and is shown to be quantized when 
the Fermi energy lies inside a gap of the spectrum. We demonstrate that the analogue of the integer quantum 
Hall effect for neutral atoms survives the non-Abelian coupling and leads to a striking fractal phase diagram. 
Moreover, this coupling induces an anomalous Hall effect as observed in graphene. 



PACS numbers: 03.75.Lm,67.85.Lm,73.43.-f 

I. INTRODUCTION 

Ultracold atoms in optical lattices offer unprecedented pos- 
sibilities of controlling quantum matter and mimicking the 
systems of condensed-matter and high-energy physics (HE). 
Particularly fascinating is the possibility to study ultracold 
atoms under the influence of strong artificial Abelian and 
non-Abelian "magnetic" fields. The experimental realization 
of artificial Abelian "magnetic" fields, which reproduce the 
physics of electrons in strong magnetic fields, is currently 
achieved through diverse schemes: for atoms in a trap the 
simplest way is to rotate the trap |2, 3], while for atoms in 
optical lattices this can be accomplished by combining laser- 
assisted tunneling and lattice acceleration methods ||4][5][6], 
using laser methods employing dark states 00, using two- 
photon dressing by laser fields [9, 10 1, by the means of lattice 
rotations HU [12] [13] [14), or, last but not least, by the im- 
mersion of atoms in a lattice within a rotating Bose-Einstein 
condensate (BEC) 1TT31 . Several phenomena were predicted 
to occur in these arrangements such as the Hofstadter "but- 
terfly" [16 1 and the "Escher staircase" [ 5 1 in single-particle 
spectra, vortex formation J2] [TT] [T7), quantum Hall effects 
16] [T3] [T8l QjO, as well as other quantum correlated liquids 
(201. 

As shown by one of us in Ref. ETI (for an alternative 
method employing dark states see Ref. 11221 ). it is simple 
to generalize the scheme of Jaksch and Zoller for generating 
artificial Abelian "magnetic" fields [4| in order to mimic ar- 
tificial non-Abelian "magnetic" fields. To this aim we have 
to consider atoms with more internal states ("flavors"). The 
gauge potentials that can be realized using standard manipu- 
lations, such as laser-assisted tunneling and lattice accelera- 
tion, can have practically arbitrary matrix form in the space of 
"flavors". In such non-Abelian potentials, the single-particle 
spectrum generally depicts a complex structure termed by one 
of us Hofstadter "moth" 11211 . which is characterized by nu- 
merous extremely small gaps. The model of Ref. ETI has 
stimulated further investigations, including studies of nontriv- 



ial quantum transport properties l23l . as well as studies of 
the integer quantum Hall effect (IQHE) for cold atoms [19], 
spatial patterns in optical lattices 1171 . modifications of the 
Landau levels [24|, and quantum atom optics Il25ll26l . 

One should note, however, that the (7(2) gauge potentials 
proposed in Ref. [21 ] and used in most of the following works 
are characterized by non-constant Wilson loops: atoms per- 
forming a loop around a plaquette undergo a unitary trans- 
formation which depends on one of the spatial coordinates. 
Although such gauge potentials are interesting per se, the fea- 
tures characterizing the Hofstadter "moth" result from this 
spatial dependence of the Wilson loop, rather than from their 
non-Abelian nature. Indeed, the Hofstadter "moth"-like spec- 
trum may actually be found in the standard Abelian case with 
a Wilson loop proportional to x (see Fig. 1). 

Two of us have shown that cold fermionic atoms trapped 
in optical lattices and subjected to artificial "magnetic" fields 
should exhibit an IQHE |fl9l . If a static force is applied to 
atoms, for instance by accelerating the lattice, the transverse 
Hall conductivity gives the relation between this external forc- 
ing and the transverse atomic current through the lattice. It 
has been shown that this transverse conductivity is quantized, 
cr X y = — j-, where C is an integer and h is Planck's con- 
stant. Note that this quantity can be easily measured from 
density profiles, as shown recently by Umucalilar et al. 11271 . 
The quantization of a xy occurs, however, only if the Fermi en- 
ergy of the system is located inside a gap of the single-particle 
spectrum. While the observation of the IQHE seems to be ex- 
perimentally feasible in Abelian "magnetic" fields, it is hardly 
so in the deeply non-Abelian regime in which the gaps of the 
"moth" become very small 11281 . 

The question therefore arises whether the consideration of 
non-Abelian gauge potentials characterized by a constant Wil- 
son loop could stabilize the spectral gaps and guarantee the ro- 
bustness of the IQHE in ultracold fermionic gases and whether 
an anomalous IQHE, as observed in graphene, can exist in 
such systems. 
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Figure 1: Energy spectrum _E = £■($) in the case of the Abelian 
gauge potential A — (0, 27r$m 2 , 0) with x — ma, corresponding 
to the non-constant Wilson loop W(m) = e l2 ' r *( 1 + 2m ). Compare 
with the Hofstadter "moth" depicted in Fig. 1 of Ref. [ 19 1 or in Fig. 
4 of Ref. |29|. The energy E is expressed in units of the tunneling 
amplitude t. 

In this work we provide affirmative answers to these ques- 
tions by considering the IQHE in a system that features a 
non-Abelian gauge potential characterized by specific non- 
commutating constant components and by a constant Wilson 
loop. We calculate the energy spectrum and we obtain a ro- 
bust band structure with well developed gaps, which differs 
drastically from the case of the gauge potential of Ref. |2"TI . 
In particular, we note the existence of van Hove singularities 
in the density of states and we obtain their analytical expres- 
sion. We then evaluate the conductivity a xy for neutral cur- 
rents using topological methods: we express a xy in terms of 
the topologically invariant Chern numbers associated to each 
energy band [30]. We eventually present a salient fractal phase 
diagram which represents the integer values of the transverse 
conductivity inside the infinitely many gaps of the spectrum. 
In this way, we show that the atomic transport is ruled by 
the quantum Hall effect supremacy: the IQHE survives in the 
non-Abelian regime, but undergoes strong modifications with 
striking similarity to the anomalous IQHE in graphene PD : 
the transverse conductivity suddenly changes sign due to the 
presence of van Hove singularities and is, under certain con- 
ditions, anomalous because of conical energy spectra. 

II. ARTIFICIAL GAUGE POTENTIALS IN OPTICAL 
LATTICES 

We consider a system of non-interacting two-component 
fermionic atoms trapped in a 2D optical square lattice of unit 



length a, with sites at (x = ma,y = no), with n.m inte- 
gers. The non-interacting limit can be reached using Feshbach 
resonances, or simply at low densities. The optical potential 
is strong, so that the tight-binding approximation holds. The 
Schrodinger equation for a single particle subjected to an arti- 
ficial gauge potential then reads 

m—l.n ) 

+ t b (u y Vw+i + u li> m.n— l) — E 1prn,n-> (i) 

where U x [resp. U y ] is the tunneling operator and t a [resp. 
tb] is the tunneling amplitude in the x [resp. y] direction. In 
the following, we use a as the length, and t a = tb = t as the 
energy units, and set h = c = e = 1, except otherwise stated. 
The tunneling operators are related to the gauge potential ac- 
cording to U x — e lAm [resp. U y = e lAy ]. 

Here we consider the gauge potential 

27r$m + 0a x ,6), (2) 

where a and j3 are parameters, (a x , a y ) are Pauli matrices and 
<I> is the number of Abelian magnetic flux quanta per unit cell. 

In order to realize such a potential one may consider the 
method of Ref. ETl . However, the specific form of this gauge 
potential allows us to consider an even more practical scheme 
based on a generalization of the method currently developed 
by Klein and Jaksch [ 15 1, or the one developed by the NIST 
group [9, 10 1 . A more detailed description of various propos- 
als can be found in Appendix A. 

The tunneling operators are 2 x 2 unitary matrices, 

U x — cos a + iu y sin a, 

U y (m) = e l2 ^ m (cos (3 + ia x sin (3), (3) 

which act on the two-component wave function V'm.n- 

The single-particle Hamiltonian is invariant under trans- 
lations defined by the operators T!J? ipm^n = 4>m+q,n and 
T y ipm.n = ipm.ji+i under the condition that <f> = |, where 
p and q are integers. Consequently, the system is restricted to 
a q x 1 super-cell and one can express the wave function as 
"0m, n = e lkxTn e lkyn u m , with u m a g-periodic function. The 
wave vector k belongs to the first Brillouin zone, a 2-torus 
defined as k x € [0, =y] and k y € [0, 27r]. The Schrodinger 
equation ([T} then reduces to a generalized Harper equation 
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III. NON-ABELIAN GAUGE POTENTIALS 

Artificial gauge potentials generally induce the following 
non-trivial unitary transformation for atoms hopping around a 
plaquette of the lattice: 

U=U x U y (m+l) U x U\ (to) . (5) 



If one sets a = dm or /3 = dir, where d is an integer, 
the matrix U — exp(i27r<I>) is proportional to the identity 
and the system behaves similarly to the Hofstadter model 
COD. When a = dn/2 and (3 = d'n/2, where d and d' 
are odd integers, one finds that U = — exp(i2-7r<I>) and the 
system is equivalent to the 7r-flux model in which half a flux 
quanta is added in each plaquette 1311 . In these particular 
cases where U = ± e i2jr *, our problem decouples into two 
independent Schrodinger equations describing particles in 
an Abelian magnetic field: the system is in the Abelian 
regime. For any other values of the parameters a and j3, the 
matrix U is a non-trivial U (2) matrix, namely U git. 1, and 
the system exhibits the non-Abelian Aharonov-Bohm effect: 
the system is in the non-Abelian regime. In Appendix B we 
present a more rigorous discussion of the gauge symmetry in 
our model, and explain in more details the meaning of the 
"genuine Abelian" and "genuine non-Abelian" cases. 

A criterion for distinguishing between the Abelian and the 
non-Abelian cases is obtained by introducing the gauge in- 
variant Wilson loop, 

W = tr U x U y (to + 1 ) Ul U\ (m) 

= 2e 427r * (cos 2 a + cos 2/3 sin 2 a). (7) 

The unitary matrix U, Eq. |6]), is proportional to the 2x2 iden- 
tity matrix if and only if \W\ = 2 (see Appendix B). There- 
fore the non-Abelian regime is reached under the condition 
that \W\ ^ 2. Besides we note that the Wilson loop possesses 
the symmetry W(a, (3) — W(f3, a). In Fig. |5j where we show 
the Wilson loop's magnitude as a function of the parameters, 
\W\ = \W(a,P)\, we can easily identify the regions corre- 
sponding to the Abelian (\W\ — 2) and to the non-Abelian 



(4) 
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In the presence of the gauge potential Eq. Q, atoms perform- 
ing a loop around a plaquette undergo the unitary transforma- 
tion: 



(6) 

I 

regimes (\W\ ^ 2). In this figure, the parameters evolve in 
the range a, (3 S [0, tt], and we note that the Abelian 7r-flux 
regime is reached at a singular point, a — (3 = tt/2. 

We also point out that the criterion according to which the 
non-Abelian regime is reached when [U x , U v ] ^ 0, and which 
can be found in previous works |fl9l |2TI . doesn't apply in 
the present context: for the situation where a — dn/2 and 
(3 = d'ir/2, where d and d' are odd integers, one finds that 
[U x , U y ] = 2i e 2l " i7r<l 'a z , while the system is Abelian because 
of its trivial Wilson loop, \W\ = | - 2e 2l7r *| = 2. 

Contrary to the non-Abelian systems considered in previous 
works lfT9ll2Tll23l . we emphasize that the gauge potential Eq. 
(|2]i leads to a Wilson loop which does not depend on the spa- 
tial coordinates. In the following section, we show that this 
feature leads to energy spectra and fractal structures which 
significantly differ from the Hofstadter "moth" ETI . 



IV. THE ENERGY SPECTRUM 

The energy spectrum can be obtained through direct diago- 
nalization of Eq. (ffl. 

In the Abelian regime corresponding to a = dn or /3 = dn, 
where d E Z, one finds q doubly-degenerated bands for 
$ = | . In this particular case, the representation of the spec- 
trum as a function of the flux $ leads to the fractal Hofstadter 
"butterfly" [16|. For the other Abelian case a = j3 = \, the 
system behaves according to the 7r-flux lattice: the spectrum 
E = E(Q) depicts a Hofstadter "butterfly" which is contained 
between $ = [0.5; 1.5], i. e. shifted by $ = 0.5 with re- 
spect to the original "butterfly", and the system remarkably 
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Figure 2: (Color online) Wilson loop's magnitude as a function of 
the parameters | W \ — \W(a, /3)|. The Abelian regime is determined 
by the criterion \ W\ = 2: in the range a,f3 £ [0, 7r], the system is 
equivalent to the Abelian Hofstadter model along the lines a = 0, n 
or (5 — 0, 7r, and is equivalent to the Abelian 7r-flux model at the 
singular point a = /3 = n/2. For any other values of the parameters 
a, P g [0, 7r], \ W\ 2 and the system is non- Abelian in the sense 
that the unitary loop operator U is not proportional to the identity 
matrix. 



describes zero-mass Dirac particles [31 1. 

In the non-Abelian regime, which is reached for arbitrary 
values of the parameters (a,j3), the degeneracy of certain 
bands is lifted and large gaps remain, as illustrated in Fig. [3] 
For these general situations, the representation of the spectrum 
as a function of the flux $ leads to new interesting features. 
As in the Abelian case, one observes repetitions of similar 
structures at various scales. However, new patterns arise in 
the non-Abelian case, as illustrated in Fig. |4]for a — (3 — j 
and in Fig. [5] for a = 1 and (3 = 2. It is worth noticing 
that for arbitrary values of the parameters (a,/3), the spec- 
tra show well-developed gaps contrasting with the Hofstadter 
"moth" which appears in the non-Abelian system proposed in 
Ref. ETi . We further notice that the spectrum is periodic with 
period T$ = 1 and is symmetric with respect to E = and 
$ = 0.5. 

In the non-Abelian regime close to a, (3 = it/ 2, one ob- 
serves that conical intersections are preserved in the energy 
spectrum. As shown in the next section, the particles behave 
similarly to Dirac particles in this non-Abelian region and the 
system exhibits an anomalous quantum Hall effect. 

We eventually note that when the flux $ = 0, the density 
of states reveals several van Hove singularities at the energies 
E = ±2(1 + cos x) an d E = ±2(1 — cos x), where x = a , P- 
As the flux increases, these singularities evolve and generally 
merge. 




Figure 3: (Color online) Spectrum E = E(k x ,k y ) for a = 1, 
(3 = 2 and $ = 0.2. While the degeneracy of some of the bands is 
lifted, the three central bands remain doubly-degenerate. Blue inte- 
gers represent the quantized transverse conductivity ha xy for Fermi 
energies situated inside the six gaps of the spectrum. The energy E 
is expressed in units of the hopping parameter t. 




E 

Figure 4: Spectrum E = E{<&) for a = /3 = \ and $ = gfp. 
where p is an integer. When $ = 0, four van Hove singularities are 
located at E = ±(2 + y/2) and E = ±(2 - \/2). The energy is 
expressed in units of the hopping parameter t. 



V. INTEGER QUANTUM HALL EFFECT AND THE 
PHASE DIAGRAM 

We evaluate the linear response of the system described 
by Eq. Q to an external force (lattice acceleration) applied 
along the y direction and we evaluate the transverse conduc- 
tivity a xy using Kubo's formula. Following the method of 
Ref. ||l9ll . one can generalize the well-known TKNN expres- 
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sion 11321 to the present non-Abelian framework, yielding 

- {9kyU X j\d kx uxj)jdk, (8) 

where is the jth component of the wave function corre- 
sponding to the band E\ such that Hu\ — E\u\, and T 2 
refers to the first Brillouin zone of the system. The Fermi en- 
ergy Ep is supposed to lie within a gap of the spectrum. The 
transverse conductivity is then given by the contribution of all 
the states filling the bands E\ < E-p situated below this gap. 

Eq. <j8j conceals a profound topological interpretation for 
the transverse conductivity based on the fibre bundle theory 
ll33l . In the present framework, such bundles are conceived as 
the product of the parameter space T 2 with the non-Abelian 
gauge group U(2). This product space, which is supposed to 
be locally trivial but is generally expected to twist globally, is 
characterized by the non-Abelian Berry's curvature 

T = (d kx A y - d ky A x + [A x ,A y })dk x dk y , (9) 

where (-4^)^- = (u\i\dk ti u\j) is the Berry's connection. The 
triviality of the fibre bundle is measured by the Chern number 

C(E X ) = f trJF, (10) 

which is a topological invariant and is necessarily an integer. 
Note that each band E\ is associated to a specific fibre bundle, 
on which a Chern number is defined. One eventually finds that 
the Hall-like conductivity Eq. ([8]) is given by a sum of integer 
Chern numbers, 

E X <E F 

As a consequence, the transverse Hall-like conductivity of the 
system evolves by steps corresponding to integer multiples of 
the inverse of Planck's constant and is robust against small 
perturbations. 

The evaluation of these topological invariants leads to a 
complete understanding of the IQHE which takes place in the 
present context. The aim is then to compute the Chern num- 
ber associated to each band E\ of the spectrum. This com- 
putation can be achieved numerically thanks to an efficient 
method developed by Fukui et al. | 34 1 and which can be ap- 
plied to our specific system. This method is summarized as 
follows: the Brillouin zone T 2 , defined by k x e [0, ^] and 
k y e [0, 2ir], is discretized into a lattice constituted by points 
denoted k[ = (fc x ;, k y i). On the lattice one defines a curvature 
T expressed as 

^ia(fcj) = lnt/i(fc,)tf a (fci + i)C7i(fc, + 2)- 1 C/ 2 (fc i )- 1 , 

(12) 



where the principal branch of the logarithm with — it < 
J~i2/i < 7r is taken, p, is a unit vector in the direction \x, 
and 



U^ki) = £>Aj(*OI«Ai(fc| + fx))/^(ki), (13) 

i 

defines a link variable with a normalization factor A/" M (fe;) 
such that \U^(ki)\ = 1. The Chern number associated to the 
band E\ is then defined by 



1 



This method ensures the integral character of the Chern num- 
bers and holds for non-overlapping bands. In the situations 
where the spectrum reveals band crossings, a more general 
definition of the link variables U^(ki) has been proposed in 
Ref. 041 . 

We first compute the Chern numbers for the specific case 
illustrated in Fig. [3] For a = 1, = 2 and $ = 0.2, 
the Chern numbers associated to the seven bands are respec- 
tively {-1; -1; 3; -2; 3; According to Eq. (fTTj), the 
transverse conductivity's values associated to the six gaps are 
{1; 2; —1; 1; —2; —1} as shown in Fig. [5]and in Fig. |6]a. 

The phase diagram describing the IQHE for our model can 
eventually be drawn. In this diagram we represent the quan- 
tized transverse conductivity as a function of the Fermi en- 
ergy E-p and flux Here we illustrate a representative exam- 
ple of such a phase diagram which was obtained for a = 1, 
(3 = 2 (cf. Fig. B). This striking figure differs radically from 
the phase diagrams obtained by Osadchy and Avron in the 
Abelian case [35| since the Chern numbers associated to the 
gaps are no longer satisfying a simple Diophantine equation 
11361 . Consequently, the measurement of the transverse con- 
ductivity in this system should show a specific sequence of 
robust plateaus, heralding a new type of quantum Hall effect. 
In order to give an alternative representation of this effect, 
we show two sections through this phase diagram, <j xy = 
a xy (E F ), at $ = 0.2 and $ = 0.02 (see Fig. ^. 
This new effect is comparable to the IQHE observed in Si- 
MOSFET or the anomalous IQHE observed in graphene in 
the "low flux" regime $ 1 corresponding to experimen- 
tally available magnetic fields. In this regime, the quantized 
conductivity evolves monotonically by steps of one between 
sudden changes of sign across the aforementioned van Hove 
singularities (see Fig. [5]and Fig. |6]b). Moreover, in the vicin- 
ity of a, (3 = 7r/2, the quantized conductivity increases by 
double integers because of Dirac points in the energy spec- 
trum, in close similarity with the anomalous IQHE observed 
in graphene. 
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Figure 5: (Color online) Spectrum E = E(<&) and phase diagram for a = l,/3 = 2 and $ = | with g < 97. Cold [resp. warm] colors 
correspond to positive [resp. negative] values of the quantized conductivity. Purple corresponds to a null transverse conductivity. For <E> <C 1, 
the quantized conductivity evolves monotonically but suddenly changes sign around the van Hove singularities located at E ~ ±1 (see the 
alternation of cold and warm colors). The Fermi energy is expressed in units of the hopping parameter t and the transverse conductivity is 
expressed in units ofl/h. 



The main experimental challenge consists here in combin- 
ing several established methods into one experiment: laser as- 
sisted tunneling |4|, BEC immersion 1 15 1, and density profile 
VI. CONCLUSIONS measurements G3. 



Summarizing, we have proposed in this paper how to real- 
ize in cold atomic systems a textbook example of non-Abelian 
gauge potential characterized by a constant Wilson loop. Our 
main result is that despite the coupling between the different 
"flavor" components of the single-particle wave functions, the 
spectrum exhibits well-developed gaps of order of 0. 1-0. 2t, 
i.e. about 50-100 nK. 

The IQHE survives in the deeply non-Abelian regime and 
acquires a unique character specific to the non-Abelian nature 
of the gauge fields. It is characterized by a particular sequence 
of robust plateaus corresponding to the quantized values of 
the transverse conductivity. Moreover, the non-Abelian cou- 
pling induces controllable van Hove singularities as well as 
an anomalous Hall effect, similar to the effect induced by the 
hexagonal geometry in graphene. Experimental observation 
of this distinctive effect requires to achieve T smaller than the 
gaps, i.e. of order of 10-50 nK, which is demanding but not 
impossible. 
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Figure 6: (Color online) Sections through the phase diagram illus- 
trated in Fig. [5] o xy = a xy (E F ), at $ = 0.2 (a) and $ = 0.02 
(b). To each plateau a xy — constant, is associated a colored rect- 
angle based on the line a xy = 0. The color code and the parame- 
ters (a = 1, P = 2) are the same as in Fig. [5] When $ = 0.2 
(a), one observes the sequence already represented in Fig. 3: the 
transverse conductivity associated to the six gaps is respectively 
{1; 2; -1; 1; -2; -1}. When $ = 0.02 (b), the quantized conduc- 
tivity evolves monotonically but suddenly changes sign around the 
van Hove singularities located at E ~ ±1. The Fermi energy is 
expressed in units of the hopping parameter t and the transverse con- 
ductivity is expressed in units of 1/h. 



A. EXPERIMENTAL REALIZATIONS 

The present work is devoted to the study of single-particle 
physics of cold atoms subjected to artificial gauge potentials. 
Recently several proposals have been made in order to real- 
ize such fields through diverse experimental techniques. For 
completeness we present in this appendix several feasible pro- 
posals for realizing Abelian and non-Abelian gauge potentials 
for cold atomic gases captured in a single trap or in an optical 
lattice. 

a. Abelian gauge fields Here we list the most promis- 
ing proposals for generating Abelian gauge potentials: 

• Rotating traps Perhaps the simplest and most common 
way is to use rotating traps. Under rotation of the sys- 



tem, the centrifugal forces act as a Lorenz force for the 
atoms dynamics. In harmonic traps, when the rotation 
frequency approaches the trap frequency, the system ex- 
pands and the atoms behave as an electronic system 
subjected to a constant magnetic field directed along 
the rotation axis |fl] El EJ. This method, which has 
been employed by many experimental groups, suffers 
from instabilities as the rotation frequency approaches 
the value of the trap frequency. 

• Rotating lattices In order to simulate an artificial mag- 
netic field for particles evolving on a lattice, one has 
to superimpose a corotating lattice on a rotating Bose- 
Einstein condensate (BEC) 11141 . In such arragements, 
the quantum Hall effect and the physics of vortices can 
be subtlely explored ifTTI fT2l fT3l fT4ll . 

• Laser methods employing dark states Perhaps among 
the most promising, these methods can be used in sin- 
gle traps and in optical lattices. This method is based 
on the fact that the adiabatic motion of a A-type 3 -level 
atom creates a non-degenerate dark state in the presence 
of external laser fields. If such a dark state is space de- 
pendent, the atomic wave function aquires a topological 
(or Berry's) phase, which can be interpreted as the ef- 
fect of an "artificial" magnetic field QUI. 

• Laser assisted tunneling and lattice tilting Several pro- 
posals have been made in order to realize "artificial" 
magnetic fields in optical lattices by using laser assisted 
tunneling |4] |5] . The first method along thoses lines 
was proposed by Jaksch and Zoller [4 1 and can be sum- 
marized as follows. The method uses two superim- 
posed lattices: one lattice traps atoms in one internal 
state (red), while the other traps atoms in another in- 
ternal state (blue). Consequently, the entire lattice is 
realized in such a way that it consists of arrays of alter- 
nating columns (red-blue-red-blue-...) along the y di- 
rection. One can use for this aim the states from differ- 
ent hyperfine manifolds. Then by tilting the lattice in 
the x-direction, one creates a linear energy level shift. 
The tunneling along the y-direction is standard, while 
the tunneling along the x-direction in laser assisted - 
it occurs via a stimulated two-photon Raman transition 
tuned resonantly to the energy shift (constant for neigh- 
boring columns). The Raman laser has a global phase 
that depends on the y-coordinate, and which induces a 
phase to the tunneling matrix elements. Again, these 
phases modify the atoms dynamics and can be inter- 
preted as the effect of an "artificial" magnetic field. 

• BEC immersion Very recently, a novel and elegant pro- 
posal has been made by Jaksch and Klein [15|. The 
method consists of immersing the optical lattice and the 
atoms of interest (^4-atoms) into a BEC consisting of 
B-atoms. The A-atoms collide with B-atoms and cre- 
ate phonon excitations, which in turn will induce effec- 
tive interactions between A-atoms [15|. If the immer- 
sion occurs into a rotating BEC, the resulting effective 
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interactions induce space-dependent phases to the tun- 
neling matrix elements of the v4-atoms. Therefore these 
phases modify the dynamics of A-atoms and simulate 
the presence of a magnetic field. 

• Two-photon dressing by laser fields Eventually one can 
use a two-photon dressing field in order to create an 
artificial vector potential, as demonstrated recently by 
Lin et al. (TO | for Bose condensed 87 Rb atoms in 
the F = 1 manifold. The created dressed states are 
momentum and spin superpositions, and the effective 
Hamiltonian corresponds exactly to the one of charged 
particles in the presence of a "magnetic" vector poten- 
tial. The magnitude of this "artifical" field is set by the 
strength and detuning of the dressing Raman fields. So 
far the method has allowed to generate a uniform vector 
potential (i.e. zero magnetic field), but the extension to 
a non-uniform potential with non-vanishing curl should 
be straightforward ifTUl . 

b. Non-Abelian gauge fields Some of the methods 
listed above and allowing for the realization of "artificial" 
Abelian gauge potentials can be extended in order to create 
non-Abelian fields. We list below the methods that can be 
considered in order to create a gauge potential of the form 
Eq.|2]), which is of interest in the present work. We focus on 
the possibilities to use alkali atoms with several internal states 
(which may be states from different and/or from the same hy- 
perfine manifold). 

• Laser methods employing dark states As pointed out in 
Ref. Il22ll . these methods can be generalized to the non- 
Abelian case, employing degenerate dark states. Such 
states arise in a four level system with three levels cou- 
pled to a common fourth (the so-called "tripod" con- 
figuration) under pairwise two-photon-resonance con- 
ditions. In particular, such dressing fields can be applied 
to atoms trapped in optical lattices. 

• Laser assisted tunneling and lattice tilting As men- 
tioned above the method of Jaksch and Zoller can be 
directly generalized to the non-Abelian case [21 1. This 
can be acheived by considering atoms in two internal 
hyperfine manifolds (red and orange) and using Zee- 
man states as "flavor" or "color" states. The case of 
U(2) gauge fields is presented in Fig. [7] The lattice 
is tilted in the x and y directions and a static magnetic 
field is used in order to split the Zeeman states energies. 
The tunneling in both horizontal and vertical directions 
is laser assisted, and the Raman lasers are tuned in such 
a way that when an atom in a given Zeeman sublevel 
tunnels, it changes its state into a coherent superposi- 
tion of all the Zeeman states. Controlling the Raman 
lasers strength and detuning allows therefore to control 
the values of the parameters a and j3 in the potential 
Eq.|2|. Additionally if the laser acting in the y direction 
has a phase that depends linearly on the x coordinate, a 
non-vanishing Abelian gauge field, characterized by the 
flux <T>, is created. 



x/a 
m + 1 



m - 1 



U y (m+J) 

3 3^3 3 a 3 

u 1 iff ^ h 

3 3^3 n 3 n 3 
I > K 3 3 



n - 2 n - 1 n n + 1 n + 2 



y/o 



V/V . 




n- 1 



n + 1 



y/a 



Figure 7: (Color online) Optical lattice setup for (7(2) gauge po- 
tentials: Dark red empty semicircles [resp. dark red filled semicir- 
cles] denote atoms in states \g\) [ resp. [32)] and light orange empty 
semicircles [resp. light orange filled semicircles] denote atoms in 
states I ei ) [ resp. |e2)]. Top: Hopping in the x direction is laser 
assisted and allows for unitary exchange of colors; it is described 
by the same unitary hopping matrix U x for both \gi) and je,:} states. 
Hopping along the y direction, which is described by the hopping 
matrix U y (m), is also laser assisted and attains m-dependent phase 
factors. Bottom: Trapping potential in the y direction. Adjacent 
sites are set off by an energy A due to the lattice acceleration, or 
a static inhomogeneous electric field. The lasers are resonant 
for transitions \gi) — > |e») for n — > rail. Because of the spa- 
tial dependence of Q±i (running waves in the ±sc direction), the 
atoms hopping around the plaquette get the unitary transformation 
U = Ul{m)U x U v {m + l)Ul 



• BEC immersion As mentioned in the text, the immer- 
sion method can be used directly to generate the poten- 
tial Eq.|2|. We consider atoms with two internals states 
(two different hyperfine states, or two Zeeman states of 
the same hyperfine manifold). The rotating BEC creates 
the Abelian magnetic field with controllable flux $ via 
phonon induced effects as demonstrated in Ref. ifTSl . 
We additionally use the laser assisted tunneling and lat- 
tice tilting methods (and if needed a magnetic field) in 
order to realize transitions that transfer the used inter- 
nal states into their superpositions. In this way we cre- 
ate non-trivial hopping operators U x ,U y <G U(2) and 
we may control the values of a and (3. In this scheme 
we note that the lasers no longer need to have a space- 
dependent phase, since the Abelian part of the gauge 
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potential <f> is induced by the rotating BEC. 

• Two-photon dressing by laser fields Finally, the method 
of dressing by Raman fields ifTol could be easily gener- 
alized in order to create a non-Abelian vector potential 
in an optical lattice. The vector potential which is lo- 
cally created should then depend on the internal states 
that are used. 

For the concrete realization of the specific gauge potential 
Eq. |2]) that we consider in this work, we may use 40 K atoms 
in F — 9/2 or F = 7/2 hyperfine manifolds, or 6 Li with 
F = 1/2. For 40 K one should optically pump and restrict the 
atomic dynamics to the two lowest Zeeman sublevels in each 
of the hyperfine manifolds. One can then employ different 
lattice tiltings in the x and y directions to perform laser (Ra- 
man assisted) tunnelings that change the internal states of the 
atoms. As mentioned above, this allows to control the param- 
eters a and j3, and fixes the non-Abelian part of the tunneling 
matrix elelments. In order to generate the Abelian part of the 
potential, i.e. the m-dependent phase factor in front of U y (m) 
in Eq. ([3]), there are two possible approaches. One possi- 
bility is to follow the lines of Ref. ET1 : the Raman lasers 
responsible for the tunneling in the y direction should have a 
x-dependent phase. The other idea is to use the immersion of 
the system in a rotating BEC 1151 . Both approaches will allow 
to control In experiments, one routinely reaches the values 
of (laser assisted, or direct) tunneling rates in the range of 5- 
10 kHz (~ 0.5 /iK), Fermi temperatures of the same order, 
and temperatures T ~ 0.2 Tp ~ 50-100 nK (see for instance 

Gill). 



B. GAUGE SYMMETRY 

In this Appendix we discuss in more details some aspects 
of the gauge symmetry proper to the models considered in 
this paper. In particular the distinction between the Abelian 
and non-Abelian cases is clarified through the introduction of 
the gauge-invariant Wilson loop. 

c. Classical external gauge fields ( CEGF) Let us start 
by making a clear statement that the gauge fields considered 
in this paper are classical, external and do not correspond to 
dynamical variables. The situation is thus analogous to that 
one encounters in atomic physics when one considers interac- 
tions of matter with electromagnetic (EM) fields (i.e. Abelian 
gauge fields), under the following assumptions: 

• The EM fields can be regarded as classical, i.e. quantum 
fluctuations can be neglected. 

• The EM fields can be regarded as external, and the back 
influence of the matter on the EM fields is negligible. 

Such situations typically occur when one deals with strong 
static electric or magnetic fields, or with strong oscillating 
fields, like the laser or maser fields. 

In the present framework, the question of interest is the study 
of quantum dynamics of matter in such fields. Such situa- 
tions are also considered in the context of non-Abelian gauge 



fields, and are widely studied in gauge field theory textbooks 
B3l . Here we are interested in the quantum dynamics of mat- 
ter (quarks) in a given (classical and external) background of 
the gauge fields. Interesting examples of such fields range 
from the fields of constant strength, to the fields carrying topo- 
logical character, such as Polyakov monopoles l43l . 

d. Gauge invariance in Abelian CEGFs: the continuum 
case We first consider the case of U(l) gauge fields (i.e. for 
instance EM fields) in the absence of the lattice and we focus 
on the non-relativistic limit. For our purpose we consider that 
the gauge potentials give rise to static magnetic fields only, 
and we assume the scalar (electrostatic) potential to be zero 
(i.e. vanishing electric field). Generalizations of the follow- 
ing concepts for more general gauge potentials are straightfor- 
ward. 

The second quantized Hamiltonian describing particles 
subjected to a U(l) CEGF reads: 

iJ=^- J dx^{x){-ihV - ^A(x)) 2 ip(x) 

+ ^ J dx^(x)^(x)^(x)^(x), (15) 

where ip(x), ip^(x) are annihilation and creation operators of 
the (single component) matter field, A(x) is the vector poten- 
tial, and g is the coupling constant for the contact interactions. 

The Hamiltonian's spectrum is invariant with respect to 
gauge transformations, in which the matter fields undergo a 
local unitary transformation, 



tp(x) 



ip(x), 



and the vector potential transforms according to 

ch. 



A(x) -> A(x) H V6(x). 

e 



(16) 



(17) 



In this context, the magnetic field B = V x A defines 
a gauge invariant field. Note that the interaction term in Eq. 
( fT5j ) is invariant under the gauge transformation Eq.(fTo*]l. Gen- 
eralizations for other kinds of gauge invariant interactions are 
also straightforward. 

e. Gauge invariance in non-Abelian CEGFs: the contin- 
uum case One can easily extend the above considerations to 
the situation in which the components of the gauge potential 
are expressed as non-commutating unitary matrices. The sec- 
ond quantized Hamiltonian describing particles subjected to a 
U(N) CEGF reads 

H=— [ dx^Hx)(-ihV - -A(x)) 2 ^(x) 
2m J ■ c 



dx ^(x)^(x)^(x)^(x), 



(18) 



where *&(x), *&i(x) are now annihilation and creation op- 
erators of the A^-component matter fields, A(x) is the non- 
Abelian vector potential, and g is the coupling constant for 
the contact interactions. The interaction term is now invariant 
with respect to the U(N) transformations of the matter fields 
which are described below. 
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Now, the Hamiltonian's spectrum is invariant with respect 
to U(N) gauge transformations, in which the matter fields un- 
dergo a local unitary U (N) transformation and the vector po- 
tential undergoes the gauge change, 



-> T(x)V(x), 



(19) 



ich r 



A(x) -> f{x)A{x)f\x) + — f{x)Vf\x). (20) 

e 

The corresponding non-Abelian "magnetic" field is defined 



by 



B(x) = curlA(a;) - —A(x) x A(x) (21) 

ch 



and transforms according to 

B(x) -> T(a;) J B(a;)f t (a;). 



(22) 



Note that the magnetic field no longer describes a gauge in- 
variant quantity in this context. 

/ Gauge invariance in Abelian CEGFs: the lattice case 
The concept of gauge invariance can also be investigated on 
lattice systems ATI l42l . The Hubbard Hamiltonian that de- 
scribes particles evolving on a lattice and subjected to a gauge 
field reads: 



" • h.r. 



(i,3) 



(23) 



where hi, b\ respectively annihilates and creates a particle on 
the site i, the phase factors exp(i(f>ij) are associated with the 
tunneling from the nearest neighboring sites j — > i, t is the 
tunneling amplitude and V is the coupling constant for the 
on-site interactions. Here we consider the case of a two- 
dimensional square lattice, but the following discussion is 
valid in general. 

In lattice gauge theories, the gauge transformation is nat- 
urally defined as an on-site U(l) transformation of the field 
operators, 



b j ->V j = e i *'b j . 



(24) 



The Hamitonian is modified under such a transformation and 
involves new tunneling phase factors given by 



Xi ~ Xj, 



(25) 
(26) 



Note that the latter expression, Eq. ( |26] l, is recognized as a 
discrete form of the gauge transformation Eq. ( fT7| i. 

As in the continuum case, the Hamiltonian's spectrum is 
left invariant under the gauge transformation Eq. (|24]l-([26]l. 
When a particle hops around an elementary plaquette, denoted 
□, it undergoes a unitary transformation \ip) — > W(D)|^)), 
where 



W(D) 



>^2(n) p i023(n) p i034(n)>j>4i(n) 



(27) 



and where the plaquette's vertices are situated at neighboring 
sites j = 1,2,3,4. This quantity, called the Wilson loop, 
is therefore associated to the plaquette □ and is defined un- 
ambiguously on the lattice by orienting each plaquette in the 
same way (say clockwise). Obviously the Wilson loops are 
invariant under the local gauge transformations, Eq. ( |26] >. 
Consequently the set of Wilson loops {W(D)} determines the 
physics completely. 

Alternatively, the physics is governed by the cumulative 
phases 

$(□) = 12 (D) + 23 (D) + cf> 3i (D) + 041 (□), (28) 

that express the Aharonov-Bohm effect: <&(□) is also a gauge 
invariant quantity and represents the magnetic flux through 
the plaquette □ in atomic units. Given a set of Wilson loops 
{W(D)} (or fluxes {$(□)}) on the lattice, there exists a 
choice of phases {</>y(D)}, such that Eq. ( |27] i is satisfied 
globally, i.e. for all loops □. In fact, there exists infinitely 
many families of local phases, {0y(D)}, that lead to the 
same physics determined by the given set of Wilson loops, 
{V^(D)}. The latter statement is always valid when open 
boundary conditions are applied to the lattice. In the case of 
closed surfaces (such as lattices with periodic boundary con- 
ditions, i.e. lattices forming a torus, or symplexes lying on the 
surface of a sphere) there is an additional physical require- 
ment: the total magnetic flux through the closed surface must 
be zero (in order to verify the Maxwell equation V ■ B = 0) 
unless the system presents magnetic monopoles. For rigorous 
proofs of these statements see for instance P4"ll45ll46ll . 

g. Gauge invariance in non-Abelian CEGFs: the lat- 
tice case The Hubbard Hamiltonian that describes particles 
evolving on a lattice and subjected to a U(N) gauge field 
reads: 



H 



b[Uijbj + h.c. 



(id) 



(29) 



where 6^, b\ are the iV-component annihilation and creation 
operators of particles on site i, JJij are the unitary operators as- 
sociated with the tunneling from the nearest neighboring sites 
j — » i, and V is the coupling constant for the on-site interac- 
tions. 

The on-site gauge transformation is now defined as a U(N) 
transformation acting on the field operators, 



b', 



T 3 b V 



(30) 



where T^T — TT^ = 1. The transformed Hamitonian in- 
volves new tunneling operators given by 



U, 



(31) 



The Hamiltonian's spectrum is left invariant under the gauge 
transformation Eq. ((30|-(|3T]l. Now when a particle hops 



11 



around an elementary plaquette it undergoes a unitary trans- 
formation — > U(D)\tp), where 

u(n) = u 12 (n)u 23 (n)u 3i (a)u 41 (a), (32) 

and where the plaquette's vertices are situated at neighbor- 
ing sites j = 1,2,3,4. This loop operator is associated 
to the plaquette □ and is defined unambiguously on the lat- 
tice by orienting each plaquette in the same way and also 
by choosing the first link. Physically, the loop operator U 
expresses the non-Abelian Aharonov-Bohm effect j47l |481 
which takes place in the system. This important effect ac- 
counts for the non-Abelian character of the system and is 
therefore observable when the loop operator is not reduced to 
a phase, U = e l. In the latter trivial case, the non-Abelian 
Aharonov-Bohm effect cannot take place and the system be- 
haves as an Abelian system. 

Under a gauge transformation, the loop operator is modified 
according to 



(33) 



and therefore this operator doesn't define a gauge invariant 
quantity. One can then consider the trace of the loop operator, 

W{D) = trU 12 (n)U 23 (D)U 34 (D)U 41 (D), (34) 

which defines a gauge invariant quantity called the Wilson 
loop. The set of Wilson loops therefore determines 

the physics of the non-Abelian system. 

Given a set of Wilson loops {W(D)} on the lattice, there 
exists a choice of tunneling operators {CZy(p)}, such that 
Eq. ( |34| ) is satisfied globally, i.e. for all loops □. More- 
over there exists infinitely many families of tunneling opera- 
tors {Uij ;(□)} that lead to the same physics determined by the 
given set of Wilson loops, {V^(D)}. 

h. Genuine non-Abelian field configurations It is now 
easy to understand which field configuration are genuinely 
non-Abelian, and which are not. Let us consider the case of a 
two-dimensional lattice subjected to a U (2) gauge potential. 

In the case in which all the loop operators are trivial and 
given by complex phase factors, £/(□) = e tS ( a 'l, we can 
gauge out all the unitary tunneling operators E/y, and replace 
them by simple phase factors. The system then reduces to 
an Abelian system. For simplicity, let us consider the case 
in which the loop operators are all equal to the same phase, 
namely U(0) = e' l0 l for all the loops □. We start with the 
first plaquette and consider the local gauge transformation in- 
duced by the following unitary operators: 



u 12 


= f\ij 12 f 2 


= i, 


yrnew 
U 23 


= fl(l 23 f 3 


= i, 


jjnew 
U 34 


= flU^fi 


= i, 


u 41 


= TlUuti 


= U{U) = 



The resulting new tunneling operators are trivialized accord- 
ing to 



(39) 
(40) 
(41) 
(42) 



This trivialization procedure can be extended to the plaquettes 
that have a common bond with the first one, and so on. Note 
that if one plaquette CL is characterized by a non-trivial loop 
operator, U (□-,) ^ e l0 \, such a gauging out is impossible! 

The system studied in this paper features a constant loop 
operator £/(□) = U, which is explicitely given by Eq. (J51. 
From the above considerations, it follows that when 
U = ±e l2,r *l, we can gauge out the unitary tunneling 
operators and reduce our problem to two independent 
Schrodinger equations for the two "flavor" components. If 
U = e l27r *l, each "flavor" component evolves in a square 
lattice subjected to an Abelian magnetic field, with flux <f>. 



When U 



1, the system describes particles evolving 



on the 7r-flux lattice and subjected to an Abelian magnetic 
field, with flux $ + n. 

On the contrary, when the loop operator is non-trivial, 
[/(□) ^ e te l, no such gauging out is possible and the system 
cannot be expressed as two independent Abelian models. 
Moreover, in this case the system exhibits true non-Abelian 
effects, such as the aforementioned non-Abelian Aharonov- 
Bohm effect. The system proposed in this work can therefore 
be used for non-Abelian interferometry, in which the final ef- 
fect of an obstacle (a unitary operator acting on the particle's 
wave function) depends not only on the nature of the ob- 
stacle but also on the particle path and location of the obstacle. 

In this work, the Abelian and non-Abelian cases are distin- 
guished through a very simple criterion illustrated in Fig. [2] 
if the Wilson loop's magnitude equals two, \ W\ = |tr U\ = 2, 
then the loop operator is trivial, namely U = er\. In this 
case, we have shown that the system reduces to two indepen- 
dent Abelian systems. In order to apply this criterion, we need 
to demonstrate the following lemma: 

Lemma: 

Given a 2 x 2 unitary matrix U € U(2), then 



\trU\ 



(43) 



where tr denotes the trace and 1 is the 2x2 unit matrix. 



Proof: 



T{ = 1, 
ft = Uu 



T} = U 12 U 23 U 3i . 



(35) 
(36) 
(37) 
(38) 



If U = e w l, then of course |tr U\ = 



»2| = 2. 



We now show that if |tri7| = 2, then U must have the 



trivial form U 



jo 



1. Since U € U{2), its eigenvalues 
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Ai = e 101 and A 2 = e 102 lie on the unit circle and we 
have that 

|tr[/| = \\xT^UT\ = |Ai + A 2 | = \e 101 + e 102 1, (44) 
= y (cos 0i + cos0 2 ) 2 + (sin 0i + sin0 2 ) 2 , (45) 
= ^2(l+cos(0i- 2 )), (46) 
where the unitary matrix T e [7(2) diagonalizes the 



matrix [7. Therefore the hypothesis |trf7| = 2 is satis- 
fied if 

0i - 2 = 2™, (47) 

where n is an integer. Under this condition, T^UT = 
diag(Ai) = e l0 \, where e = Ai = A 2 . Eventually 
this shows that U = e l0 l, since TT^ = T^T = i. □ 
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